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Abstract 

We study whether a 125 GeV standard model-like Higgs boson can be accommodated within the 
scale-invariant NMSSM in a way that is natural in all respects, i.e., not only is the stop mass and 
hence its loop contribution to Higgs mass of natural size, but we do not allow significant tuning 
of NMSSM parameters as well. We pursue as much as possible an analytic approach which gives 
clear insights on various ways to accommodate such a Higgs mass, while conducting complementary 
numerical analyses. We consider both scenarios with singlet-like state being heavier and lighter 
than SM-like Higgs. With A-terms being small, we find for the NMSSM to be perturbative up to 
GUT scale, it is not possible to get 125 GeV Higgs mass, which is true even if we tune parameters 
of NMSSM. If we allow some of the couplings to become non-perturbative below the GUT scale, 
then the non-tuned option implies that the singlet self-coupling, k, is larger than the singlet-Higgs 
coupling, A, which itself is order 1. This leads to a Landau pole for these couplings close to 
the weak scale, in particular below ~ 10 4 TeV. In both the perturbative and non-perturbative 
NMSSM, allowing large A x , A K gives "more room" to accommodate a 125 GeV Higgs, but a tuning 
of these A-terms may be needed. In our analysis we also conduct a careful study of the constraints 
on the parameter space from requiring global stability of the desired vacuum fitting a 125 GeV 
Higgs, which is complementary to existing literature. In particular, as the singlet-higgs coupling A 
increases, vacuum stability becomes more serious of an issue. 
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1 Introduction 

Very recently, the ATLAS and CMS collaborations at the large hadron collider (LHC) have 
discovered a resonance with a mass of around 125 GeV which is consistent with a standard model 
(SM) Higgs boson [1], Here, we assume it is a SM-like Higgs. 

In the SM, the quartic coupling for the Higgs and hence its mass is a free parameter so that 
the above observation can be easily "accommodated" . However, the SM is plagued by the Planck- 
weak hierarchy problem for which weak-scale supersymmetry (SUSY) is perhaps the most studied 
solution. In turn, in the minimal supersymmetric SM (MSSM), the Higgs tree- level quartic coupling 
is given by the electroweak (EW) gauge coupling so that there is actually a prediction for the Higgs 
mass. The flip side is that the Higgs mass then has a tree-level upper bound of mz, which is well 
below the above observation of 125 GeV. It is well-known that the contribution of the superpartner 
of top quark (stop) at the loop-level can raise the Higgs mass in the MSSM, but reaching 125 GeV 
this way requires stops to be as heavy as ~ 5 — 10 TeV or tuned, large At [2, 3, 4]. The heaviness 
of the stop results in reintroduction of fine-tuning of the EW scale. This issue is embodied in the 
the relation 

that is one of the minimization constraints of the MSSM potential. Here the soft mass is 
largely determined by the stop mass via RGE, resulting in | > 0.1m? [5]. Given the above 
stop mass, such large value of m 2 Hu needs to be canceled by other contributions at the right-hand 
side of eq. (1), in particular, either against the supersymmetric mass for Higgs doublets (//-term) 
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or the other Higgs soft mass ) [6, 7, 8]. At any rate the need for a cancellation of the large 
rrijj signals that a new "little hierarchy" needs to be explained in the MSSM. This little hierarchy 
problem in MSSM can be quantified with a fine-tuning measure [9]: 

(flog 777-7 , x 

A = max ° z , 2 
e t dlogOi 

with A < 5, i.e., less than 20% tuning, conventionally taken as typical for a natural theory. 

The MSSM also has another drawback, partly connected with the issue of having a cancellation 
in eq. (1) involving /i. The issue is about why is the /i-term close to the weak scale as required 
phenomenologically: there is lower bound on it of ~ 100 GeV from chargino mass, whereas much 
larger values would require tuning in eq. (1). An attractive solution to this /i-problem is the scale- 
invariant NMSSM (for a review see [10]). Here, an explicit //-term is forbidden in the superpotential, 
and instead is generated by the vacuum expectation value (VEV) of an added SM gauge singlet, 
S, which is coupled to Higgs doublets, H u and Hj, via the superpotential term 

XSH u H d . (3) 

In turn, this VEV for the singlet is driven by soft SUSY-breaking mass terms for singlet. As there 
are no explicit mass scales in the superpotential, this model can be referred to as the scale-invariant 
NMSSM. This form of superpotential can be ensured by a ^-symmetry, which can be extended 
to the soft SUSY breaking terms as well. Thus it is often also referred to as the Z3 invariant 
NMSSM K 

As a bonus of the NMSSM, we get an extra tree-level quartic coupling for the Higgs doublets 
from the same singlet-Higgs coupling A which solves the /i-problem. This effect in raising the Higgs 
mass relative to its MSSM prediction is well known and has been investigated at length in the 
literature (see for instance [13, 14, 15]). Here we re-visit this issue in the light of the discovery of 
a SM-like Higgs boson at 125 GeV. 

Variations of the NMSSM may allow dimensionful terms in the superpotential. In this work, 
we focus on the version of the NMSSM with no mass scales in superpotential so as to keep its other 
merit of solving the /i-problem. We assume that the stop contribution to Higgs mass is small, as 
follows from natural mz,At. This means that we aim at using tree-level contributions to get close 
to 125 GeV, about 110 GeV or more, and use the stop contribution to gain just the "last" 15 GeV 
or less in mass 2 . 



1 The original Z3 NMSSM typically suffers from the difficulty of simultaneously solving domain wall problem and 
tadpole problem which destabilizes electroweak scale [11]. A simple way to resolve this issue is to assume a low- 
scale inflation such that the domain walls are significantly diluted [10]. Alternatively, a suitable R-symmetry can be 
imposed to constrain the form of explicit /^-breaking terms which may reconcile the tension of solving domain wall 
and tadpole problems [12]. 

2 The estimate of 110 GeV lower limit will be detailed in the following and roughly corresponds to small it — tn 
mixing and the contribution of a stop with mass of 500 GeV which is a general upper bound from naturalness 
according to eq. (2) . 
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Our goal is to see if it involves any tuning of parameters of the NMSSM to reach the required 
tree-level value for Higgs mass. Such a suspicion is rather motivated. For example, it is well known 
that in a large region of the NMSSM parameter space, the lightest CP even state, s±, is identified 
as the SM-like Higgs. The mass of si has an upper bound at tree-level 

m 2 Sl < m| cos 2 2/3 + A V sin 2 2/3 . (4) 

If we prefer to preserve the conventional unification of the gauge couplings, which is a major success 
of MSSM, A is constrained to remain perturbative up to the GUT scale. This means that the value 
of A at the weak scale which enters the above bound cannot be much larger than 0.6 — 0.7 [16]. 
In turn, even for such a maximal value of A, the tree-level upper bound for the Higgs mass is very 
close to 110 GeV. Given that Higgs mass of 110 GeV or more at tree- level is our goal, it is clear 
that we must actually saturate the upper bound in this case. However, the upper bound is not a 
generic value of the Higgs mass in the NMSSM, i.e., we expect that the model parameters must 
be arranged in a specific manner to fit the data in this case. Here we shall investigate in detail 
what are the conditions that define the locus of maximal Higgs mass in this region of the parameter 
space of the NMSSM. 

A key quantity to identify the regions that lead to the largest Higgs boson mass (when it is the 
lightest CP-even state) is the mixing of the doublet-like states with the singlet-like states. Indeed 
it has been known that the extra quartic coupling of the NMSSM for the SM-like Higgs is not the 
"end of the story" since mixing with the singlet modifies the mass of the SM-like state [13] and in 
general pushes it away from the value of the right-hand side of eq. (4). When the singlet- like scalar 
is heavier, as discussed so far, its mixing with the doublet-like state results in a "pull-down" of the 
mass of the latter. The crucial point is that, even though the singlet-like scalar is heavier than the 
SM-like Higgs, it turns out that the shift of the mass of the SM Higgs-like state due to mixing with 
singlet state can still be significant, since this effect does not quite decouple even if singlet VEV is 
large 3 . As mentioned already, if the singlet-Higgs coupling A is to remain perturbative up to GUT 
scale, then it turns out that the maximal tree-level SM-like Higgs mass - before negative mixing 
effect - can only be barely around or lower than our target of 110 GeV. Thus one has to work hard 
to minimize the negative mixing effect on the SM-like Higgs mass in order to obtain 125 GeV mass 
at loop-level. 

There is actually a different possibility, namely, that the singlet-like scalar is actually lighter 
than the SM-like Higgs. In this case, their mixing provides a "push-up" effect, which raises the 
SM-like Higgs mass. Thus, the tree- level mass, before the mixing is taken into account, can be below 
110 GeV, lessening the tension with perturbativity mentioned above for the other case. However, 
we face a different potential worry: the state lighter than 125 GeV necessarily has a component 
of Higgs doublet so that it couples to the Z boson and a bound from LEP2 applies. Thus, in this 

3 Although we might actually require (S) ~ dew in order to have a natural /i-term, \{S). 
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case, one has to work hard to make sure that the singlet-like scalar is not too light or too strongly 
coupled to the Z boson. 

As we proceed with our above analysis, we distinguish two cases: the "perturbative" NMSSM 
that we mentioned above, which requires all the couplings do not hit a Landau pole till GUT scale, 
and the non-perturbative NMSSM, where such singularity of the couplings can appear below the 
GUT scale. The motivation is clear: the former case can manifestly maintain the successful gauge 
coupling unification of the MSSM. Another useful "axis" of our exploration of the parameter space 
is the size of soft trilinear SUSY breaking terms for Higgs doublets and singlet, A K and A\. We 
consider both cases where they are small and cases where they are large. This division clearly 
is relevant for some SUSY breaking mediation schemes that predict the A-terms to be small, for 
example, minimal gauge mediation. Finally, as far as possible, we try to come up with analytical 
insights into these various issues. 

The issue of the Higgs mass in the NMSSM has been studied recently, also in connection with 
the experimental finding of a 125 GeV SM-like Higgs boson. However our study is significantly 
different in a number of aspects with respect to the current literature. For example in Ref. [3] an 
explicit mass terms for the singlet and Higgs doublets are allowed. Other studies such as [17, 18] 
are mostly numerical and they lack a clear division into the above cases that we highlight. 

Our results can be summarized as follows, beginning with the case with negligible ^4-terms. For 
the perturbative case, we demonstrate that even if we allow tuning of parameters of NMSSM it is 
not possible to reach 125 GeV mass: for the sub-case here of pull-down, this limitation is due to 
a combination of the facts that the singlet-Higgs coupling A is small to begin with and that the 
reduction in Higgs mass due to the mixing cannot be made small. On the other hand, for push-up, 
the failure to reach 125 GeV Higgs mass is due to the light singlet-like state typically being ruled out 
by LEP2. In the non-perturbative NMSSM, the singlet-Higgs coupling A and singlet self-coupling 
k can be larger. Thus in this case the maximum tree level value itself can be well above 110 GeV. 
Typically, one then uses the singlet-Higgs mixing to get down from the above bound to 110 GeV 
at tree-level. In this way we completely avoid the tuning associated with requiring the pull-down 
mixing effect to be near vanishing that plagued the perturbative case above. Looking closer at 
the non-perturbative case it turns out that, if one wishes to avoid any kind of tuning of NMSSM 
parameters, then we end up in a specific region, namely, singlet-Higgs coupling, A, is of order one 
and singlet self-coupling, k, is larger. Such a choice implies a Landau pole for these couplings which 
is rather close to weak scale. Finally, with A-terms, for both perturbative and non-perturbative 
cases, there is more freedom, for example, it is possible to minimize the size of mixing in order to 
reach 125 GeV even in the perturbative case, but in this case it still requires tuning of A\ against 
/I. 

One of the constraints in addition to rrih ~ 125 GeV that we explicitly take into account is the 
stability of the vacuum. Although this is a complicated issue to deal with, it is essential for the 
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phenomenological viability of the model: the parameter choice that accommodates the right Higgs 
mass should correspond to a cosmologically stable vacuum. This issue, however, is either overlooked 
or only investigated in a partial way in most existing literature trying to explain rrih ~ 125 GeV 
within NMSSM. In our work we try to do a more careful study in this regard, both analytically 
and numerically. 

Here is the outline of the rest of the paper. We begin with reviewing the model and giving 
relevant formulae for our analysis. This is followed by a general discussion of how to get to Higgs 
mass of 125 GeV, including relevant issues such as perturbativity and vacuum stability. After this 
setting of the stage, we discuss in detail the viability of the model in the push-up scenario first and 
then the pull-down case, before concluding. 



2 The Model 

The scale-invariant NMSSM is a very well-known extension of the MSSM which was originally 
proposed to overcome the problem and can also serve to lift the tree-level Higgs mass. Nonetheless, 
in this section, for the sake of completeness and to fix our notation, we give some details of the 
model relevant for our discussion. For a complete overview on the NMSSM we refer to [10]. 
The superpotential of the model is 

Wnmssm = XSH u H d + -S 3 , (5) 

where S is a singlet chiral superfield and H u = H®) and H d = (H%, are the MSSM Higgs 
doublets. The soft-SUSY breaking Lagrangian is 

-Aaft = m 2 H jH u \ 2 + m 2 Hd \H d \ 2 + m 2 s \S\ 2 + XA x H u H d S +^kA k S 3 . (6) 

Hence the complete scalar potential in this sector is: 

^scalar — ^soft ~i~ ^ ^ 

{\Hl? + \Ht\ 2 - \H°\ 2 - \H~\ 2 f + | \H+H°* + H°H- d *\ 2 , (7) 

where the last line is the MSSM D-term, g2 ~ 0.65 and g\ « 0.35 denote the SU(2)l and U(l)y 
gauge couplings, and F{ = dW/dcfti are the F-terms. 

We expand the neutral scalar fields around the vacuum expectation values (VEVs) as follows 
(see, for example, [10] ) : 

H° u = v u + 4= [ (h° v + iG°) sin (3 + (H° + iA° v ) cos p] , (8) 
V2 L J 

H° d = v d + ^=[{K -iG°) cos p-(H%-iA°) sin /?] , (9) 
S = s + ^={h s + iA° s ) , (10) 
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where v = \Jv% + v^ « 174 GeV, tan/3 = v u /vd and G° is the neutral would-be Nambu-Goldstone 
boson. 

The extremization conditions of the scalar potential are 

2 2 

ml u+ ^ + X\l + 9 -l±^{vl- v ^] - Vd »B = 0, 



4 

2 2 

2 i 2 i \2 2 i 5l + 92 I 2 2\ 
+ X v u + ^ («d - ««) 



u u //B = , 



s (m| + kA k s + 2k 2 s 2 + A 2 v 2 — 2n\v u Vd) — Xv u VdA\ = , (11) 

where [i = As and i? = A\ + «s. We can use the above three equations in order to replace the 
parameters , mji d and m| by the three VEVs. 

The CP-even mass squared matrix in the above basis of h®, H®, h® is 

(X 2 v 2 sin 2 2/3 + m| cos 2 2/3 2rw 2 cot2/3 2A 2 su - 2t> 2 i? \ 

-2x; 2 r + 2A ^+ 2 A / s -2Rv cot 2/3 , (12) 

ks(4ks + ,4 k ) + |jA A Asin2^ / 

where r = - fj) sin 2 2/3, and R = 4 A(ks + \A X ) sin 2/3. 

The basis used in eq. (12) has the advantage that it shows clearly what state couples linearly to 
the SM gauge bosons, such that it can be produced at LEP2 in association with a Z boson. In fact 
h® is rotated in the same way as v, so that it is exactly the linear combination of H® and H% which 
is responsible for the masses of the W and Z gauge bosons. Thus only h® has tree-level couplings 
to WW/ZZ, while H® and h° s do not. This implies that only the component of h® in each CP-even 
mass eigenstate is relevant for the LEP2 limit. 

Similarly, to the extent that ATLAS and CMS collaborations at the LHC have observed a SM- 
like Higgs at 125 GeV, in particular with SM-like coupling to WW/ZZ, this state (denoted by h) 
should be predominantly h®, with small admixture of H®, h®. 

The CP-odd scalar mass squared matrix in the above basis of A® and A® is given by 

M 2_(^^ Xv(A x -2 KS ) \ 

I • A^(A A +4K S )sin2/3 _ I ■ ^ 

The Lagrangian is such that we get a light CP-odd scalar in the two approximate symmetry limits: 
the U(1)r symmetric limit when the A-terms are small; and the Peccei-Quinn symmetric limit 
when k and A K are small. The light CP-odd scalars then correspond to pseudo-Nambu-Goldstone 
bosons from spontaneous breaking of these symmetries. If their mass is below about 10 GeV, there 
can be a constraint coming from T-decay into these CP-odd scalars. In general, the CP-odd states 
can also be produced in association with the CP-even scalars at LEP2. 
Finally, there are physical charged Higgs bosons whose mass is given by 

m H± = — + v V - A • V 14 ) 



sin 2/3 V 2 
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Although it is not directly relevant for our analysis, we refer the reader to [10] for the mass 
matrix and the properties of the neutralinos and charginos. 

Without loss of generality we take A > while n and (i can have either signs. Meanwhile, in 
the examples shown in the following sections we take k > and \i > 0, as we do not expect the 
other choices of signs to result in any qualitative change of our findings. 

3 Getting SM-like Higgs of mass 125 GeV 

As mentioned above, the 125 GeV SM-like Higgs discovered at the LHC is likely to be predominantly 
h®, with a small admixture of H®, h®. Neglecting this mixing, its mass squared is given by A / (f 1 
in eq. (12): the m| piece is the same as in the MSSM and the A 2 piece results from the extra 
quartic of the NMSSM. However, as mentioned in the introduction, the mixing of h® with H® and 
hg modifies the mass from the above value, lowering (raising) if mixes with heavier (lighter) 
states. In fact the SM-like Higgs mass can be written as: 

™l = m z cos 2 2(3 + Sm\ loop + 

A 2 i> 2 sin 2 2/3 + <5m 2 mix . (15) 

We further discuss below, in turn, each term in the above formula. The first line is the MSSM 
value, including the stop mixing effect for which we use the approximate formula from eqs. (69-71) 
of [19]: 

Mr 



2 _ 3m f 

0m h loop — 2tt 2 v 2 



'»S + T + ^( 3m?/ " 2 - 16 ^( X ' + 21 "S)] (16) 



where the "stop-mixing" parameter Xt is given by 

2 (A t - n/tirnp) 2 r [A t -ju/tan/3) 5 



X t 



M 



•i L l2M i 



(17) 



fht = mt/ (1 + 4a s / (3tt)) is the MS top quark mass, m t being the pole mass, g s is the QCD coupling 
and M~ = m^m^ is the geometric average of the two stop mass eigenvalues. 

We do not assume ad-hoc large mixing in the stop sector which tends to increase the Higgs mass. 
We also assume 500 GeV as a rough upper limit of stop mass Mf as it is favored by naturalness of 
weak scale . So we take typical values /i = 200 GeV, tan/3 = 1.5, < 500 GeV and At = in 
eq. (16) and find the stop loop contribution to rrih can be up to 15 GeV. Therefore our goal here is 
to achieve a tree-level mass for the SM-like Higgs of 110 GeV or more, up to 125 GeV. Of course 
if one would allow a heavier, hence less natural stop our analysis would be overly restrictive. 



4 As emphasized in reference [3], for A ~ 2, even 1.5 TeV stop mass can be natural. However, for such large values 
of A, the tree-level bound is already well above 125 GeV and so a large loop contribution to Higgs mass from a heavy 
stop is not really needed anyway. 
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Moving onto the second line in eq. (15), the first term here is the contribution characteristic of 
the NMSSM coupling A. Notice that for a moderate or large A, in the NMSSM, nih is maximized 
at tan/3 ~ 1 in contrast to the case in the MSSM, where nih is maximized at tan/3 — > oo. Thus, 
we choose small tan/3 in our analysis. As it will play an important role in what follows, we have 
introduced the notation <5m? mix in the second line here to encode the effect from the mixing of h® 
with H% and h° s . 

For the discussion here, and for our analytic study in Sections 4 and 5, we shall neglect the 
h® — H® mixing and only deal with the h® — h° s mixing part of Sm\ mix , reducing the 3-by-3 mixing 
problem into a simpler two states mixing problem. The justifications of this approximation are as 
follows. Typically we have to choose parameters so that the h® — H® and h® — H® mixings are 
small effects. The reason is that H® must be somewhat heavy, because H® is accompanied by a 
charged Higgs of similar mass which is constrained by b — > s'y [23]. As H® gets heavier its effect 
on rrih decouples since the h® — H® mass mixing term does not scale with H® mass, as shown in 
eq. (12). Furthermore in the NMSSM small tan/3 is well motivated and the h® — H® mixing actually 
vanishes at tan/3 = 1. In contrast, the effect of the h® — h® mixing on SM-like Higgs mass cannot 
be decoupled. The main reason is that both the h° v - h° s mass term and h® mass scale with the 
VEV for S. 

Similar to the h° - H% mixing we neglect the hP a — H® mixing since it affects the mass of the 
doublet like state only at higher order and it also vanishes at tan/3 = 1. We stress that we neglect 
h® — H® and h® — H® mixing only to give analytical arguments, of course, the full mixing is included 
in our numerical analysis. 

As anticipated in the introduction, there are two possibilities for the mass of the SM-like Higgs 
boson as far as the effect of the mixing effect in eq. (15) is concerned: 

• The "pull-down" case where the singlet-like is heavier than SM-like Higgs, such that the 
mixing reduces SM-like Higgs mass. 

In this case, (i) M\ x < -M§ 3 and (ii) > (110 GeV) 2 such that after the pull-down mixing 
effect, we end up with a SM-like Higgs tree-level mass of 110 GeV. Obviously, eq. (12) then 
suggests that A > 0.7. 

• The "push-up" case where the singlet is lighter that the SM-like state, so that the mixing 
increases the SM-like Higgs mass. 

By definition this implies that (i) M\ x > M% 3 and (ii) M\i < (110 GeV) 2 since, after the 
effect of the mixing is included, we want to get to 110 GeV at tree- level. Therefore, based on 
eq. (12), we must have A < 0.7. 

Elaborating on the above discussion we are in position to outline the following argument. If the 
model has to be perturbative up to the GUT scale, then in the pull-down scenario, the region of 
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Figure 1: The green solid line is the maximal mass of the SM-like Higgs boson for tan/3 = 1.5 from 
eq. (15) in the "pull-down" scenario. The two horizontal blue dashed lines are the reference values 
125 GeV, that is about the observed mass of the SM-like boson, and 110 GeV, that is the minimal 
value of the tree-level mass to attain 125 GeV after the inclusion of radiative corrections from the 
stop, taking M 4 ~ < 500 GeV. The light (dark) purple shaded regions mark the values of A that, 
according to the one-loop RGE, get larger than v2~7T, our definition for non-perturbative, below the 
GUT scale for k = 0.1 (k = 0) at the weak scale. To allow larger values of A to be compatible with 
perturbative unification we have introduced extra matter in 4 copies of 5 + 5 of SU(5) at 5 TeV. 

the NMSSM parameter space that can give a SM-like Higgs mass as large as 125 GeV is not at all 
a generic one. For the push-up case we shall indicate how we get a similar conclusion, though for 
a different reason. 

If the coupling A must remain perturbative up to the GUT scale, it cannot exceed a generously 
estimated upper-bound at about 0.7 . Even plugging in such large value of A in eq. (15) one 
finds that, in the pull-down scenario, effect of the mixing must be minimal to keep the physical 
Higgs mass above 110 GeV at the tree-level. This fact was already mentioned in the introduction, 
but now we demonstrate it in Figure 1 and it suggests that the parameters of the model must be 
arranged in a specific pattern such to have 

6m 2 hmix ~0. (18) 

This conspiracy of parameters seems rather unnatural, or can be taken as a suggestion for a peculiar 
relation that must come from a UV construction, and therefore we introduce a measure of "tuning" 
to quantify to what extent we are requiring cancellations among apparently unrelated parameters 
of the model. 

In a fashion similar to what is done for the tuning of the EW scale (or Z boson mass) as given 
in eq. (2), we introduce the tuning measure for the Higgs mass and expand the expression using 

5 In fact A is typically constrained to be below or around 0.6, as in general the need for a non- vanishing k reduces 
the maximal value of A compatible with perturbativity up to the GUT scale. 
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eq. (15) and (12): 



A. 



max 

6i 




■mix 



a! log Oi 



max 

6i 



dlog {ml -M^) 



max 

Oi 



d\og6i 

d\og (m? — X 2 v 2 sin 2 2/3 — m?7 cos : 



2 2/3) 



(19) 



d\og9i 



where 0j = {A, «, A\, A k , fj,} are the various parameters of the NMSSM evaluated at the weak scale, 
where we compute the Higgs mass. This measure of tuning tries to capture how tightly related 
must be the parameters of the NMSSM to attain 125 GeV for the mass of the SM-like Higgs state. 
A large value of A m j x corresponds to a fine tuning among the values of two or more parameters. 

In the following sections we will discuss both analytically and numerically what relations among 
the NMSSM parameters are implied by the recent LHC observation of a SM-like Higgs with mass 
125 GeV. Furthermore we will quantify using eq. (19) how natural are the identified regions of the 
parameter space. Altogether we will give the coordinates of the islands in the natural region of the 
NMSSM parameter space spotted by the LHC Higgs result. 

For the push-up scenario the discussion follows a rather different path, but comes to a similar 
conclusion. The push-up scenario has a singlet-like state at the bottom of the spectrum. This 
state can be within the reach of LEP2 and to keep it unobservable in the dedicated searches of 
LEP2 experiments it must be significantly decoupled from the Z boson, i.e., it must be mostly a 
pure singlet. In conflict with this is the fact that the lift of the Higgs mass, Smf imix , goes to zero 
as the mixing angle vanishes. Consequently in the push-up scenario the need for an almost pure 
singlet-like state at the bottom of the spectrum and the need to lift the Higgs mass are in tension. 
This picks out special regions of the parameter space where both the LHC and LEP2 constraints 
can be accommodated. In a spirit similar to that of the case of pull-down we can quantify how 
"special" is the choice of the NMSSM parameters to attain such decoupled singlet in association 
with the needed lifting of the tree-level mass of the SM-like state. Since the issue is also about 
obtaining a small 5m 2 imix , we can use the same measure of tuning introduced for the pull-down 
case as given in eq. (19). 

In the following sections we will organize our analytical and numerical discussion on how to 
attain the mass of the SM-like Higgs at 125 GeV according to the push-up versus pull-down clas- 
sification discussed above, paying particular attention to the size of A m i x . Before getting into this 
detailed analysis, we would like to comment on two other issues of the NMSSM, which become 
especially relevant for such a large mass of the SM-like Higgs. Namely, the stability of the desired 
vacuum and the presence of Landau poles for the couplings, i.e., perturbativity of the model up to 
high scales. 
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3.1 Perturbativity 

As we noted when commenting about eq. (15), tan/3 must be "small" in order to make use of the 
extra quartic to increase the SM-like Higgs mass beyond its MSSM value. In fact for moderate 
or large A as used in this case, the tree-level upper bound on Higgs mass is itself maximized at 
tan/3 ~ 1. However, such low tan/3 implies that the top Yukawa coupling is larger than 1 already 
at the weak scale. In addition to the danger that top Yukawa hits Landau pole below GUT scale, 
such a choice, in turn, makes A blow up faster 6 . Therefore we shall choose tan/3 = 1.5 as a 
"representative" value for our analysis. 

Also, in order to ameliorate the potential problem with Landau poles below GUT scales, we allow 
the possibility for extra matter at low scales. Such matter should appear in complete representations 
of SU(5) so as to preserve gauge coupling unification, and tends to give larger MSSM gauge couplings 
in the UV. The gauge couplings, in turn, prevent top Yukawa and singlet-Higgs couplings from 
blowing up too quickly. 

3.2 Constraints from unrealistic minima 

A set of important constraints on the values of the NMSSM parameters is that the ones fitting 
nih ~ 125 GeV should correspond to a viable realistic vacuum. Basic checks include requiring 
it to be a minimum with all scalar masses squared being positive. On top of this the desired 
vacuum there should be either a global minimum or a meta-stable minimum at cosmological time 
scale. Because of more parameters in NMSSM scalar sector, the check of vacuum stability is more 
sophisticated than in MSSM, in general with no neat analytic condition. Existing works often 
overlook the issue of the global stability that, in particular, may become a crucial one when A is 
large [24, 25]. In this work we take this issue seriously and conduct a thorough analysis on how 
vacuum stability may constrain the parameter space that fits ~ 125 GeV. In this section we 
first give a general analytic formulation for such an analysis. Then, we show a complete analytic 
study for a particular yet practically interesting limit where A4f 3 = 0, which, in turn, leads to 
^ rn h mix ~ 0- R- ecau that this is the condition which was identified as the key to attaining the 
desired Higgs mass while preserve perturbativity up to GUT scale: see the discussion after eq. (15). 

Guided by this analytic understanding, in the numerical results contained in the following 
Sections 4 and 5 we will impose vacuum stability constraints on the parameter space that has an 
otherwise acceptable vacuum, i.e., which gives ~ 125 GeV and, is compatible with the LEP 
bound when it is relevant. In the numerical studies the assessment of the global or local nature of 
the desired minimum is carried on by a global and unconstrained numerical search of possible extra 
minima of the potential. This must be contrasted with the checks performed by standard tools [26] 
that evaluate analytical expressions for a (in general incomplete) set of the extremal points of the 

6 On the other hand negative stop 2-loop contribution is a bit smaller in size due to a cancellation between the 
now larger top Yukawa coupling and QCD coupling, see eq. (16). 
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one-loop effective potential. Although limited by the precision of the numerical scan and by the 
fact that we use the tree-level potential, our method is in a sense more general as it can capture 
extra minima in any direction in field space. As we will see later, the requirement of absolute 
vacuum stability may rule out significant portions of the parameter space that are allowed by all 
other constraints. Our result gives a warning regarding the viability of some regions of parameters 
space and motivates a more detailed analysis at loop-level, that we leave for future work. 

We start from stationary conditions of the NMSSM potential which determine the VEVs at 
extremum points, as given in eqs. (11). In general these 3 equations are coupled and there are 6 
extremum solutions: 

• ( s = 0,u = 0), (s = 0,v = v U0 ); 

• (s = s T , v = 0), (s = s r , v = v r ); 

• (s = s Ul ,v = 0), (s = s Ul ,v = v Ul ), 

where the "r" label denotes the desired quantities at the realistic vacuum with rrih ~ 125 GeV, 
v 2 = m 2 Zr /g 2 , g 2 = (gf + g^)j1 and s T = fi T /X corresponding to natural /x r . On the other hand, 
the "u" label denotes the nonzero quantities in extra unrealistic vacua. With A / there is no 
runaway direction in NMSSM potential. Therefore in order to decide whether the desired vacuum 
is a global minimum, it is sufficient to check if V T is deeper than all other 5 extremum points. As 
already mentioned, in order to keep our analytic study at a manageable level, we will work with 
the tree-level potential. We are aware that, due to the importance of loop effects, this may not be 
accurate enough in some situations. In our numerical results presented later we shall consider the 
comparison of minima as "unresolved" when the value of the potentials at the two points differ by 
less than 5% 7 . 

Applying stationary conditions on H u ,Hd one obtains two practically interesting variables 
sin 2/3 e g , m 2 z eS for each of the above solutions as functions of input parameters and the corre- 
sponding fj, eS : 

Sm2f3eS = m 2 + m 2 + 2u 2 +AV ' ( } 

m H u + m H d + ^cff + A V eS 



ml, eS = == . Hd -m 2 Hu -m 2 Hd -2M 2 . (21) 

(sm2/3 eff ) 2 



Here in order to clarify things, we use the label "eff" as a general notation for output parameters 
from stationary conditions. Input parameters are those intrinsic to a model: A, k, A\, A k , rrVjj , 



7 This threshold of 5% is our rough estimate of size of EW loop effects, which are expected to play a significant role 
in breaking the tree-level degeneracy of the two vacua. Some numerical investigation of the impact of loop-corrections 
computed using NMSSMTools [26] confirm that our estimate is in the right ballpark. In fact the smallness of loop 
correction is expected once one chooses a natural stop mass. 
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m H d > m s- Of course one set of (sin 2f5 e s , m z eff , ^ e ff) needs to be the values at the desired, realistic 
vacuum, which as mentioned we will label with the label "r" . 

Fixing X, k, A\, A K we can retrieve the mass squared model parameters from the quantities at 
realistic vacuum 

" sTn2ft 2Mr XV * (22) 



'2(A A + f Mr ) Mr 2 
+ I — ^ m Z ,r 



sin 2/3, 



r 



A\ 2 ^ Vl- (sin 2ft) 



*»k - 2(A : i :^ > "' - 2 ^- A ' (23) 

m| = ( A 2 s r u 2 + 2k 2 s^ — Xkv 2 s t sin 2/3 r — -AA\t> 2 sin 2/3 r + k^4 k s 2 

Then one can (in principle) solve for other 5 extremum solutions (s, v) and their corresponding 
(sin 2/3 e ff , 

m z eff' Meff)- Finally we can plug all the relevant quantities into the potential and compare 
the depths at different extrema. To achieve this goal we rewrite the potential as: 



4 2 m z,es( sin 2 M m| eff (cos 2/3 eff ) 
16g 4 



V**n = " A 2 ,,.„, +V> m , (2 1) 

where 



Klin — ^Meff + 2 ^j^K/^cff + T2 m lMeff • (25) 

Here a few comments are in order. Based on the first two terms in eq. (24), among solutions with 
the same /i e g, the one with nonzero mz,eS always has deeper potential. Therefore, among the 5 other 
extremum solutions we mentioned before, it is actually sufficient to just check if (s = 0, v = v UQ ) or 
(s = s Ul , v = v m ) has a deeper V than desired vacuum. Therefore later in this section we focus on 
these two unrealistic minima. 

The solution of (s = 0, v = v no ) is a special case, which is easily obtained from eqs. (21) and (20): 

sin2/3| s=0 = 0, (26) 
m z \ s=0 = -2m% u , 

where we make the conventional ansatz rn? H ^ < m^. This choice is motivated by radiative sym- 
metry breaking, however choosing rrijj^ > n r ijj d does not change the result. 

On the other hand, the solution (s = s Ul ,v = v Ul ) is hard to solve because this requires 
fully solving the 3 coupled cubic equations. However, there is a special limit where doublets can 
decouple from the stationary point equation of S, which greatly simplifies the process of obtaining 
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the solution. To see this we recast the stationary point equation of 5, the last one of eq. (11), in 
the form 



2kV + A K s + m 2 s s + Xv 2 (Xs - us sin 2/3 A\ sin 2(3) = 0. 



(27) 



Apparently, doublet VEVs can be eliminated from eq. (27) when v = or the expression in the 
parenthesis vanishes. Remarkably this requirement coincides with .M13 = which in turn implies 
the small singlet-doublet mixing eq. (18). As we discussed before this is necessary to attain a 
125 GeV mass for the SM-like Higgs while being compatible with perturbative unification. 

In our following discussions we will first consider the competition between realistic vacuum and 



the generic (s = 0, v 



solution given in eq. (27), in particular, the dependence on model 



parameters. Then we conduct a complete analysis for the simple special case mentioned above, i.e., 
when doublets decouple from stationary equation of S. 



3.2.1 V r vs. V u (s = 0,v 



At the realistic vacuum, ^ e g- = fj, T , mz : eS = m z,r, sin2/3 e g 
expression of mg in eq. (24) into eq. (24), and simplify V r as: 



sin2/3 r . One can then plug the 



2 m Z. 



r (sin2/3 r 



m^ r (cos2/3 r 



1% 4 8g 2 
2/i r - sm2(3 r (A x + ^Mr) 



A 4 



3A 3 



(28) 



Notice that the second line in eq. (28) vanishes in the limit when the doublets decouple from the 
stationary point equation for S. 

At the unrealistic extremum point (s = 0, v = v UQ ), we clearly have V^ in = in eq. (24). Thus, 
the only relevant quantity here is mz, e s which can be written in terms of quantities at the realistic 
minimum ju r , m^,, sin2/3 r using eqs. (24,27). Finally, we get the following potential at this point: 



V u (s = 0,v = v uo ) 



m 



Z,cff 



(29) 



. ^±M^ + 2 Mr 2 + A\ 2 



+ 



sin2/3 r 



2{A X + 
sin 2/3 r 



\ 2 v 2 Vl- (sin 2ft) 2 



So, the question of whether the realistic minimum is the global one or not reduces to a compar- 
ison of the value of the potential in eq. (28) and eq. (29). We first make some general qualitative 
comments that will be later confirmed by our quantitative studies. With all the other parameters 
fixed, at a value of ^ much larger than mz, the unrealistic extremum tends to have deeper V than 
V r . A general "hint" of this feature, that may also apply to the other unrealistic extrema, comes 
from the fact that the rn^ eS term(s) in eq. (24) give negative contribution to V. At the realistic 
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vacuum mz, e ft has the value 91 GeV, achieving which, according to eq. (21), requires tuning when 
fj, r is large. On the other hand, at the other extrema with non-zero m^j, the size of mz t eS is 
typically mz,u ~ m H u ~ f-h- Therefore V u gets a stronger negative push from larger | — m A z eff | 
terms. To be more specific to this extremum (s = 0, v = v Uo ) we can see from eq. (28) and eq. (29) 
that, at large fi r it is the dominant negative contribution of 0(fj>f) in the potential which tends to 
make V u deeper. In fact V n D — ^ [l — f (1 — cot 2/3 r )] 2 ^ to be confronted with the contribution 
in the potential at the realistic minimum V T D — fiA'r- With tan/3 r ~ 1.5, A ~ 0.7 as a relevant 
example, we find the coefficient of /x 4 term in V u is larger in magnitude (but negative) than in V T , 
except for a small intermediate region of 0.5 < k < 1. 

Furthermore one can study the dependence on A with all other parameters fixed. At larger A, the 
unrealistic extremum tends to have deeper V than V r . This can be seen by comparing the leading A 

Ttx^ (sin 2/3 )^ 

terms in the two potentials. These are z ' r 16g4 r A 2 C V r and — (g^(l — yl — (sin2/3 r ) 2 ) 2 t> 4 A 4 + 

cfi 2 v^X 2 ) C Vu where the c— term has an O(l) coefficient depending on Clearly, the leading 

0(A 4 ) term gives a significant negative contribution to V u which is absent in V T . 

We can also study the dependence on k with all other parameters fixed. At larger k, the 
unrealistic extremum tends to be shallower than V r . This can be seen by comparing coefficients of 
the leading k dependent term of O(k) at large k. These are — ^/J*f C V T and — - c V u . 

With tan f3 T ~ 1.5, A ~ 0.7 as an example, V r gets a larger negative contribution at larger k. 

3.2.2 A complete vacuum analysis at doublet-singlet decoupling limit 

As mentioned earlier, it is much easier to solve for all extrema in the NMSSM potential in the 
limit of decoupling doublet from singlet. Strikingly, this is also the interesting limit for obtaining 
rrih — 125 GeV while being compatible with perturbative unification and LEP limits. Here we shall 
show that there are parameter regions where the singlet-doublet decoupling holds and the realistic 
vacuum is a global minimum. Despite the special nature of the limit where our finding applies, 
this constitutes a proof of existence of a vacuum that is both stable and able to accommodate 
rrih ~ 125GeV. The results may be generalized to other cases as well. 

Notice that in order for our approach to be self-consistent, other extrema, i.e., including 
(s = s ui ,t> = u ui ), also needs to satisfy the decoupling condition A4 2 3 = 0, just like the realis- 
tic minimum: 



2k 

2/j, eS - sin 2(3 eS (A\ + — /x e ff) 



. (30) 



This more restrictive condition of vanishing mixing entry also disfavors tachyons at other extrema 
which makes them more "competitive" compared to realistic vacuum. Meanwhile such additional 
condition constrains input parameters A\,A K to satisfy certain relations with X,k and output 
quantities at realistic vacuum. 
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In this limit the stationary point equation of S eq. (27) reduces to 



s(2kV + A K s + m|) = , 



(31) 



and the solutions {0, s T , s 



u 



} take a simple form. Evidently s n satisfies: 



s u + s r = -A k /(2k) . 



(32) 



We solve for nonzero s u using eqs. (20,21,30). First, A\ is fixed by the decoupling relation eq. (30) 
and fi T , sin 2/3 r , A, n. Then based on eq. (21) and eq. (30) we can rewrite v e g and yu e g as function 
of sin2/3 e ff, then plug in eq. (20). A consistency equation for sin2/3 e g- only can be obtained. We 
solve this equation and pick out 3 real solutions: 0, sin2/3 r , sin2/3 u . The corresponding solutions for 
(s e g, v e s) can be easily derived. 

In order to compute the potentials we need to retrieve the remaining model parameters A K , rrig. 
We plug s u into eq. (32) to find the model parameter A K . Plugging in m^, sin 2/3 r , s T at realistic 
vacuum and A, k, A\, A k into eq. (24), we retrieve the last model parameter m|. Now we have 
obtained all 5 other sets of stationary solutions (s,v) and can compute potentials using eq. (24) 
and see when the realistic vacuum has the deepest potential, i.e., it is the global minimum. 

The result of our analytic study is shown in Figure 2, where we evaluate the value of the potential 
at all the stationary points as a function of fi r , i.e., with all the values of the input parameters 
A\, A K , m? s fixed as functions of fi T , tan/3, k, A as explained above. The red line corresponds to the 
value of the potential at the realistic realistic minimum. For different choices of k in the two panels 
this figure illustrates for what range of \x x the realistic minimum is the global one. 

In the figure we observe the danger of taking large [i r discussed in Section 3.2.1 above. The left 
panel shows that the value of the potential at the extremal point (s = 0,v = v uo ) gets deeper than 
the realistic minimum when /i r is large. The trend as we increase k is as expected from the earlier 
discussion: the realistic minimum is shallower at smaller k, and larger k is found to disfavor the 
unrealistic extremal point (s = 0,v = v Uo ). 

4 Push-up Scenario 

As mentioned, one possibility of alleviating the tension between perturbativity up to GUT scale 
and rrih « 125 GeV is the "push-up" scenario. As discussed earlier, H® is more decoupled from 
and hP s , in particular at small tan/3 as favored in the NMSSM to obtain large enough Higgs mass. 
Therefore we focus on the 2-by-2 mass matrix in the basis of (h®, h®). In this region of parameter 
space the SM Higgs-like state (mostly h®) is typically the second lightest state of the spectrum 8 , 
heavier than the singlet-like state (mostly h®), thus the mixing with the lightest state induces an 

8 We find that in the alternative case where the SM Higgs-like state is the heaviest state of the spectrum, the 
charged Higgs would be too light to be consistent with direct searches and flavor constraints such as from b — > sj. 
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A=0.7, k=0.2, tan/3=1.5 A=0.7, «:=0.8, tan/J=1.5 




100 200 300 400 500 100 200 300 400 500 



//(GeV) /j(GeV) 

Figure 2: The value of the potential evaluated at all possible stationary points defined by the 
requirement of having \i r = /x at the realistic minimum. The parameters of the potential are fixed 
as functions of /i r as discussed in the text. In red and thick line we show the realistic minimum 
(s = s r ,v = v t ). The other extremal points of the potential are detonated by dashed line with the 
following color code: orange for (s = s r ,v = 0); blue for (s = 0, v = v UQ ); brown for the origin 
(s = 0,v = 0); cyan for (s = s Ul ,v = v Ul ); green for (s = s Ul , v = 0). The orange line is close to red 
line in particular at large k as a result of the fact that the potentials of the corresponding extrema 
only differs by mzeS dependent terms in the potential eq. (24) which is a sub-leading contribution. 



increase in its mass. Thus, we must have 

M 2 n >M 2 33 . (33) 

Apparently, it seems plausible to avoid the tuning associated with a small Sm^ mix in this case. 
However, as we shall demonstrate in this section, this scenario is strongly constrained by the LEP 
bound. In the end, the allowed push-up effect is very limited which corresponds to fine-tuning 
worse than approximately 20% 9 . 

In detail, we have the following sub- mass matrix for tan/3 ~ 1: 




This tree-level mass matrix has two eigenvalues, mi and 771,2, that are the masses of the doublet-like 
and singlet-like state, respectively. In this case mi ~ m^. Also, in the push-up scenario mi > m%. 

Loop corrections to the masses from the stop sector are known to be substantial. As mentioned 
earlier, for a representative analytic study we choose m^ ~ 500 GeV and negligible At. Such mass of 
the stop gives a substantial correction to the Higgs mass yet still being compatible with naturalness. 
Including the loop correction from such stop, we can aim at the tree-level Higgs mass 

mi « 110 GeV. (35) 

9 See Ref. [27] for a similar analysis. 
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In the push-up scenario the mass of the doublet-like Higgs is increased by mixing effects, there- 
fore the range of the coupling A of interest is given by 

Xv < mi « 110 GeV (36) 

and thus 

A<0.7. (37) 

Also, to achieve the required spectrum with a lighter singlet-like state, we require the condition 
of eq. (33) which explicitly reads as 

v 2 

AV > Ak 2 s 2 + A k ks + — A X X. (38) 

The lightest CP-even Higgs of the push-up scenario is typically in the kinematical reach of LEP 
since m% < 110 GeV. This state is in general a mixture of singlet and doublet, due to the latter 
component it is constrained by LEP limits reported in Figure 2 of Ref. [28]. As we will see, a 
significant portion of the parameter space of the push-up region is excluded by LEP. 

To get an analytic understanding, we approximate the upper limit from LEP on the singlet- 
doublet mixing from above reference as a step- function of the lighter scalar mass: 

Region-I sin 2 9 w 0.01, < m 2 < 80 GeV , (39) 
Region-II sin 2 9 sa 0.1, 80 GeV < m 2 < 100 GeV , 
Region-Ill sin 2 9 « 0.4, 100 GeV < m 2 < 110 GeV , 

where 9 is the mixing angle between the singlet and doublet in the matrix eq. (34). 

Note that (for all values of 9) we have tan 29 = 2M\ 3 /{M\ l - M 2 3 ). Thus, the bound on 
mixing angle can be phrased directly as a constraint on the matrix elements in eq. (34). In turn, 
the constraint on mass matrix elements can be mapped onto constraints on the model parameters 
such as A, k etc. We can then get an idea of the tuning level using the same parametrization 
discussed in Section 2. 

The quantification of this tuning is the main goal of this section. For this purpose we shall 
compute numerically A m ; x as defined in eq. (19). However in order to facilitate the analytic un- 
derstanding, it is convenient to define a simplified fine-tuning formula in terms of mass matrix 
elements. We start from the reduced 2-by-2 mass matrix eq. (34). In case of Region-I and II 
we may expand mass eigenvalues w.r.t. small mixing angle. Furthermore to obtain analytically 
tractable expressions we need to assume M. 2 33 -C M.\i- The latter assumption is not necessarily true 
in all the parameter space, however it can be realized in large fraction of the relevant parameter 
space since the push-up condition requires M\i > ■M% 3 - With s « v 10 we get the approximate 

10 From chargino bound on ft term, we know s > v, so taking s ~ u gives a conservative estimate of the level of 
tuning here. Furthermore, as we will see from numerical analysis section, after including other constraints such as 
no-tachyons, fi term has to be < 120GeV which means s ~ v. 
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expression derived from eq. (19): 



l 13 

For later convenience we also rewrite this fine-tuning as 



A ■ ~ ^ (40) 
mix — M 2 ' 



M\ x - m\ ' 



A ■ ~ - (d'\) 



It is easy to check that this closely reproduces the A m ; x that we would obtain from the standard 
definition eq. (19). 

We first consider the small mixing region including Region I and II in eq. (40) where sin 2 6 < 0.1 
and m 2 < 100 GeV. We then have tan 2(9 < 1 so that 9 « M^/iM^- Ml 3 ). We can expand the 
mass eigenvalues in a series in 0: 

m\ = M 2 11 + 6 2 (M 2 U -Mh), (42) 
m 2 = M 2 3 - 9 2 (M 2 U -M 2 3 ). (43) 

From eq. (42) we get 

A&= (1 + (44, 
and plugging this back into eq. (43) we obtain 

o 1 + 26< 2 , 9 1 + e 2 9 

ml = Q2 M{ x ~ —Q^ m l ■ ( 45 ) 

From this we get bounds on M\i and on the degree of tuning at the boundary (i.e., maximum 
value of 9) of different mass regions 11 : 

• Region-I: sin 2 9 = 0.01, < m 2 < 80 GeV, m x = 110 GeV 
Plugging this in eq. (45) we get 



109.5 2 GeV 2 < M 2 n < 109.7 2 GeV 2 , 
24 2 GeV 2 < M 2 13 < 34 2 GeV 2 , 

ll 2 GeV 2 < M 2 33 < 80 2 GeV 2 , (46) 
that corresponds to at best about 2% tuning (i.e., A m i x w 40) from eq. (40). 

• Region-II: sin 2 9 = 0.1, 80 GeV < m 2 < 100 GeV, mi = 110 GeV 
Plugging in eq. (45) we get 

11 Tuning is clearly worse for smaller values of 9. 
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108 2 GeV 2 < M 2 U < 109 2 GeV 2 , 
24 2 GeV 2 < Mj 3 < 39 2 GeV 2 , 

83 2 GeV 2 < M 2 33 < 101 2 GeV 2 . (47) 

This would correspond to a tuning of about 5% (A m ; x = 20). However we notice that the 
obtained ranges for the matrix elements are incompatible with the assumptions that we made 
to get an analytically manageable expression for the tuning. Therefore this estimate of the 
tuning in Region-II is expected to be very crude. In fact the numerical results given later 
show that the tuning can be as good as 20% (A m i x = 5). 

Here a few comments are in order. From the above analysis can only lie in the narrow 
region 108 2 — HO 2 which directly tells A ~ 0.6, and not smaller. Additionally we remark that the 
closeness of M\\ to m\ and consequently the tuning is a result of the constraint from LEP. 

These bounds on the mass matrix entries, covering regions I and II can be translated into 
constraints on the model parameters A, k, A k , A\. The discussion is divided into several sub-cases 
depending on which ^4-terms are vanishing. This division is general from a bottom-up point of view 
and also can be easily mapped onto complete models of SUSY breaking where a typical size of the 
A-terms is predicted. 

• Case-I: A K = 0,A X = 

The small value(s) of the bound on Ai\ 3 in eq. (46) and eq. (47) combined with the expression 
for it in eq. (34) and A ~ 0.6 as derived above implies 

A ps k ss 0.6. 

We know from chargino LEP bound that fi = As > 105 GeV, thus s > v. Meanwhile 
M 2 3 = 4k 2 s 2 , hence A ss k implies -M 2 3 > M.\\ contradicting both the push-up condition 
M.\i > -A^| 3 given in eq. (38) and the A^| 3 bound that we just derived as in eq. (46) or 
eq. (47). In conclusion this region does not work in any way. 

• Case-II: A K ^ 0, A\ = 
In this case we also need 

A ss k w 0.6 

since A R does not enter M\ 3 - The push-up condition and A4 33 bound can now be accom- 
modated by a suitable choice of A K . However, we remark that this case is not viable if we 
require perturbativity up to GUT scale because k is too large. 
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• Case-Ill: A K = 0,A x ^O 

In this case the small M\, is attained requiring (A — k)s — A\/2 ~ 0. When k is taken small 
the model may be perturbative up to the GUT scale. In this case a tuning between \i and 
A\ is required. Later in our numerical studies we will find this tuning to be 20% or worse 
(A m i x > 5). The chosen A\ also needs to satisfy the bound on Ml, as derived in eq. (46) or 
eq. (47) as well as the push-up condition M\ x > A4§ 3 . 



Of course we would expect the combination of Case-II and Case-Ill (A K / 0, A\ ^ 0) may 
contain viable parameter space as well. But we leave this to our numerical study as the 
analytic study gets complicated in this combined case. 



For Region-Ill the mixing is large, tan 29 > 1 (sin 2 9 > 0.15). In this case the proper way to 
get simplified mass eigenvalues is to expand in a series in e = (M.\i — Ml,)/{2M 2 ,) < 1: 



nil 



■Mil + M 2 „ + 2-^13 U + 



4(.Mp 2 



'13; 



(l + -)^ 2 1 + (l- e --)^I 3 



Mi 1 + Mi 3 -2M( 3 [l + 



(M* u - Ml,? 



4(-M 13 ) 2 



(l--)_M 2 1 + (l + e + -).M§3 



Eq. (48) gives 



M 



2 _ 2mj-(l + \)Ml 1 
1 - e 



1 



(48) 



(49) 



(50) 



In this region, expanding the eigenvalues of 2-by-2 mass matrix w.r.t e and further simplify with 
M\i > .M33, we derived an approximate tuning formula from eq. (19): 



A, 



2M 2 U 



(51) 



Plugging this into eq. (49) and applying e = 1/5 (sin 2 9 sa 0.4), and 100 GeV < rri2 < 110 GeV, 
we get the following bounds on A^f l5 A^ 2 3 : 

• Region-Ill: sin 2 9 w 0.4, 100 GeV < m 2 < 110 GeV, m 1 = 110 GeV 



107 2 GeV 2 < M 2 U < 111 2 GeV 2 , 
34 2 GeV 2 < M 2 n < 46 2 GeV 2 , 

103 2 GeV 2 < Ml, < 109 2 GeV 2 , (52) 
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that corresponds to a tuning at best about 10% (i.e., A m i x « 10) from eq. (51). Notice 
that now Mfi and -M§ 3 have to be near degenerate, which might be seen as another kind of 
tuning. 

Some comments on the above large mixing case are in order. Despite of the large mixing, the 
bound on tells that we still need 



and thus the push-up effect is limited. This is the case because the large mixing mostly can only 
be due to the near degeneracy of the masses of the two eigenstates in order to avoid LEP limits. 
Consequently, although the fine-tuning defined by A m ; x is slightly better than the small mixing case 
of Region-I, it is still at the level of 10% or worse. As far as the discussion of tuning is concerned, 
it is worth stressing that, while the larger mixing slightly alleviate the fine-tuning of <5"i| m i x of 
eq. (15) at the same time it requires another kind of tuning such that M.\\ and M.\^ are near 
degenerate as mentioned earlier. A simple estimator of the precision of the cancellation required 
here may be (M.\i — ■A4| 3 )/(A^f 1 + AI33) - which is found to be about 5% tuning or worse. 

It is also handy to obtain constraints on other model parameters in this large mixing case. 
When the A-terms are zero, to get the mixing angle one needs 0.55 < k < 0.6, which is still very 
close to A and violates the push-up condition eq. (38), analogous to the case of small mixing, and 
thus is not viable at all. For non vanishing A-terms the discussion also mostly follows the one 
presented for the small mixing case of Region-I and Region-II. 

The analytic considerations exposed above are powerful estimators of the degree of tuning 
necessary to obtain a viable spectrum with the lightest CP-even state being sufficiently singlet-like 
to avoid LEP limits and a heavier state with a mass of 125 GeV, i.e., the SM-like Higgs boson 
observed at the LHC. To verify our findings and to cover the intermediate cases that do not clearly 
belong to one of the limits studied above, we next perform a numerical study of the properties of 
the NMSSM in the region of the parameter space that includes the cases outlined above. 

The study is performed fixing the parameters 

fj,, A x , A k , tan/3, 

while we let free the parameters A and k. All the observables are computed with no approximations. 

Following the discussion on the parametric behavior of the Higgs masses and mixing presented 
above we discuss the following points in the parameter space, in order to illustrate the cases with 
A\ 7^ 0, Ak <C A\, with A\ <C A K , A k ^ 0, and with generic A-terms, respectively: 

• n = 120 GeV, A x = 250 GeV, A k = , tan p = 1.5 

• n = 110 GeV, A\ = 50 GeV, A k = -250 GeV, tan = 1.5 
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• fj, = 120 GeV, A x = 250 GeV, A k = -250 GeV, tan^ = 1.5 

We shall identify the region in the (A , k) plane that is compatible with a number of requirements 
of the scalar potential. In particular we shall require that the the point (H u , Hj) = (0, 0) in the 
fields space is unstable, in order to favor EW symmetry breaking. For the same reason we shall 
require the potential to have a minimum at the point where the singlet field scalar field is (S) = 
s = j and the doublet fields are at the point ({Hj) , {H u )) = (v sin (3 , v cos /3), where v ~ 174 GeV. 
At the realistic minimum we must also make sure that the £7(1) of the electromagnetism is not 
broken. All in all we require the following basic conditions for a realistic potential: 



In Figure 3 we show in bright yellow, light yellow and orange regions where all these conditions 
are satisfied. The difference between these three regions lies in the additional requirements on the 
stability of the realistic minimum. Indeed, the potential of the NMSSM can have several minima 
and, as we discussed in Section 3.2, the realistic minimum may not be the global minimum, posing 
a question on the stability of the vacuum where we are supposed to live. The points in the orange 
region are points at which the potential has minima deeper than the realistic one, i.e., the realistic 
minimum is a local but not a global minimum. In the region in bright yellow the realistic minimum 
is the global minimum of the potential while in the light yellow region the realistic minimum is 
almost degenerate with another minimum. In particular in the light yellow region the two minima 
(regardless of which is the deeper) are degenerate up to a 5% difference in their depths, i.e., 



where V\2 are the depth of the two minima. This 5% threshold is chosen as we expect that loop 
corrections will change the depths of the minima by a few percent and therefore it is only when 
eq. (54) does not hold that we can use the tree-level potential to reliably determine what minimum 
is the global one. 

To firmly exclude or allow some of these points, a proper analysis would be required, namely, the 

computation of the lifetime of the vacua to establish if we are living in a vacuum that is sufficiently 

long-lived 12 . The lifetime can be estimated using the Euclidian action of the field configuration 

that interpolates between the two minima, a task that we leave for future work. However it is clear 

that the points in the orange region are likely to have too a short lifetime for the realistic vacuum 

and are in some sense disfavored. 

12 Even that would be just a necessary check for the viability of these points, since other pathologies may arise. 
For instance, there is still the risk that in the early Universe the scalar fields stay at the unrealistic minimum for too 
long and significantly alter the many successful predictions of standard cosmology. 



(H u , Hd) = (0, 0) is an unstable point 
V(v,s) < 
m 2 Si > 0, m% 1 > 0, m 2 H± > . 



(53) 



\Vi/V 2 - 1| < 0.05 



(54) 
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Figure 3: Constraints and fine tuning in the (A, k) plane in the push-up scenario. The solid 
magenta lines are the isolines of the fine tuning measure A m j x defined in eq. (19). The cyan region 
is the region compatible with both LEP searches for the lightest CP-even and with the tree-level 
mass of the SM-like state in the range 110 - 125 GeV as suggested by LHC data. The green region 
would be compatible with LHC but is ruled out by LEP. 

The light and bright yellow and the orange areas denote the regions where the realistic minimum 
is an acceptable vacuum according to the requirements of eq. (53). In the orange region other 
minima where v ^ v r or s ^ s r (or both) are significantly (more than 5%) deeper than the realistic 
minimum. In the light yellow region the extra-minima are degenerate within 5% with the realistic 
minimum. In the bright yellow region the extra minima are either absent or significantly shallower 
than the realistic minimum. The purple dashed line is the boundary of the region that contains 
couplings that stay perturbative (less than \/2n) up to Aqut = 1-6 • 10 16 GeV. 
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In Figure 3 we show in cyan the region allowed by LEP and LHC data on the SM-like Higgs 
state and singlet-like state in the push-up scenario. We also show the isolines of the fine-tuning 
computed with the exact formula eq. (19). The purple dashed line is the boundary of the region that 
contains couplings that stay perturbative (less than \/2vr) up to Kqut = 1-6 • 10 16 GeV according 
to a one-loop RGE with boundary conditions at 1 TeV and 4 sets of 5 + 5 of SU(5) entering in the 
RGE at 5 TeV. 

Figure 3 shows that, according to our estimates above, in the push-up scenario it is possible to 
have a lightest singlet-like CP-even scalar compatible with LEP and a heavier SM Higgs-like state 
at the mass suggested by the recent discovery of the LHC experiments. The cases where A\ is 
chosen to obtain a small M.\% are the least tuned, i.e., there are phenomenologically viable points 
of this type having A m j x > 5. Furthermore for sizable yl-terms it becomes possible for the model 
to be perturbative up to the GUT scale. 

Of course a detailed investigation of the second lightest CP-even scalar of the push-up scenario 
should be carried on to fully assess the compatibility of this region of parameter space with all the 
available information about the newly discovered SM-like Higgs boson. We just remark that this 
state tends to be SM-like, as it cannot mix much with the singlet, and the other state doublet-like 
state is rather decoupled. Hence the 125 GeV boson in the push-up scenario seems in reasonable 
agreement with current observations of the properties of the new boson. 

5 Pull-down Scenario 

In the pull-down scenario the effect of the mixing in eq. (15) is to reduce the mass of the SM-like 
Higgs boson relative to the upper bound. As suggested by Figure 1, the observed mass of the Higgs 
then requires A > 0.7. We divide the range of A into A ~ 0.7 and A > 0.7. The reason to do this 
splitting is that the former value can be - but is not guaranteed to be - perturbative up to the 
GUT scale, whereas the latter values are clearly non-perturbative. 

We start with an analytic study in order to have a clear view before the numerics. For this 
purpose, we assume a hierarchical mass squared matrix, i.e., M.\\ -C -^2233 an d -M12 <C 
-A/ff 3 <C A4§ 3 . Note that the above assumption are compatible with what is often referred to as 
"decoupling limit", i.e., s>u with A ~ k, as well as the case where all the states are mixed, s ~ v 
with k>A, 

As we have explained already we consider only the two state mixing problem for h® — H®. Hence 
the mass of the doublet-like state is 



where p = k/X. The first line above - which is M\i from eq. (12) - is the usually quoted upper 




m| cos 2 2/3 + A V sin 2 2/3 




(55) 
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Figure 4: Constraints and fine tuning in the (A, k) plane in the pull-down scenario. The solid 
pink lines are the isolines of the fine tuning measure A m ; x defined in eq. (19). The green region is 
the region where the tree-level mass of the SM-like state in the range 110 - 125 GeV as suggested 
by LHC data and by Naturalness considerations on the mass of the stop. 

The light and bright yellow and the orange areas denote the regions where the desired minimum 
is an acceptable vacuum according to the requirements of eq. (53). In the orange region other 
minima where v ^ v r or s ^ s r (or both) are significantly (more than 5%) deeper than the realistic 
minimum. In the light yellow region the extra-minima are degenerate within 5% with the realistic 
minimum. In the bright yellow region the extra minima are either absent or significantly shallower 
than the realistic minimum. The purple dashed line is the boundary of the region that contains 
couplings that stay perturbative (less than y/2n) up to Agut = 1.6 • 10 16 GeV. The orange dashed 
lines in the right plot are the isolines of the Landau pole scale Alp where one of the couplings gets 
non-perturbative. 



bound on SM-like Higgs mass in the NMSSM. The second line is the negative effect from mixing. 

The above formula can be generalized to include A-terms, assuming that the A-terms do not 
dominate over singlet VEVs: 



m l 



m| cos z 2/3 + A V sin" 2/3 



2„,2 



-X z v 



sm2/3(l + ^)-± 



1 _|_ Al _|_ A A sin 2/3 v 2 



(56) 



5.1 A -0.7 

In this case, M.\\ is close to (110 GeV) 2 . Thus, we cannot afford a sizable pull-down effect so that 
the goal here is to minimize this effect and saturate the upper bound on SM-like Higgs mass. This 
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minimization of the mixing is clearly a source of tuning, a tuning of the mixing as we dub it in 
order to distinguish it from other kind of tunings to be encountered later and is quantified by A m ; x 
in eq. (19). 

For negligible ^4-terms we see from eq. (55), that the only option to minimize the pull-down 
effect is to tune k vs. A such that 

/csa A/ sin 2/3. (57) 

In this case, the value of k is so large that the couplings become non-perturbative below the GUT 
scale. 

For sizable A-terms, we can keep k small so that the theory is perturbative up to the GUT 
scale. Instead, we can tune A\ vs. /i in order to minimize the effect of mixing in eq. (56): 

A\ ~ 2/i/ sin 2/3 . (58) 

Another possibility would have been to make A K large and positive so as to suppress the mixing 
effect. However for typical needed to reduce the mixing effect the CP-odd pseudo-scalar masses 
from eq. (13) are tachyonic, hence this option is not viable. 

5.2 A>0.7 

In this case the theory becomes non-perturbative below the GUT scale due to the large value of A at 
the EW scale, regardless of the value of k. The maximal Higgs mass attainable is now larger than 
125 GeV, thus a certain amount of pull-down effect is actually needed. For this reason the tuning 
of the mixing encountered for A ~ 0.7 is absent for larger A. therefore we begin the discussion with 
examples for A ~ 1. 

For negligible A-terms it is better to re-write the Higgs mass formula as 

m\ ~ "m\ cos 2 2/3 + 

AV f — =■ + - sin 2/3 | . (59) 

V p p ) 

Given the large value of A, to obtain the desired mass of the Higgs we need the term in parenthesis 
to be much smaller than 1 versus being ~ 1 for A < 0.7. For this purpose we need to deviate from 
eq. (57). There are two options. The first one is to choose p < 1/ sin 2/3 that is 

k< A/ sin 2/3. (60) 

The small k helps to keep the Landau pole far from the EW scale (even though still below GUT 
scale). However, adjusting the parameters in this way tends to require a cancellation between the 
two terms inside parenthesis in second line of eq. (59). For example, A = 1, tan/3 = 1.5 requires 
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k — 0.65 to get Higgs mass down from upper bound of 164 GeV to 125 GeV which comes from 
a cancellation of about 1 part in 5 between two terms in second line of eq. (59). Note that this 
tuning is different from the tuning of the mixing considered earlier for the case with A ~ 0.7. 
The second way to get the Higgs mass from eq. (59) is to choose p > 1, that is, 

k> A/ sin 2/3. (61) 

For this choice of parameters the tuning associated with cancellation in the second line of eq. (59) 
is not an issue since the second term always dominates. However in this case k is larger than A and 
thus perturbativity is lost rather close to the EW scale. For example, A = 1, tan/3 = 1.5 requires 
K — 3.2. For examples of a possible UV completion of the NMSSM with choice of couplings A > 1 
and k > A see Refs. [21, 22]. 

Finally, in the range of A in between 0.7 and 1 there might be natural choices of parameters. 
For example, tan/3 = 1.5, A = 0.77, k = 0.67 or 1.1, can get Higgs mass down from upper bound 
of 129 GeV to 125 GeV which is associated with a cancellation of about 1 part in 5. 

Guided by the analytical intuition developed above we now proceed to the numerical computa- 
tion of several properties of the NMSSM. As we did for the push-up scenario, in order to study the 
regions identified above we shall vary the parameters A and k while fixing the other parameters. 
The right panel of Figure 4 illustrates the case with negligible A-terms. Color codes for the global 
minimum analysis are the same as for the study of the push-up scenario in Section 4. The green 
band in the plot is the region where the Higgs mass is in the range 110 - 125 GeV. As expected 
from our analytic study, in this case it is impossible to avoid a Landau pole below GUT scale. 

Indeed the Landau pole is rather low, as illustrated by the dashed lines that correspond to 
a Landau pole at 10 4 or 10 2 TeV. If one imagines an embedding of the NMSSM in a model 
of supersymmetry breaking such as gauge mediation, where there is a meta-stable NLSP, there 
are interesting consequences from the low values of the Landau pole that we found. Indeed if 
one assumes no perturbativity breakdown below the SUSY breaking scale, then the Landau pole 
provides an upper limit on yF. Interestingly the region that accommodates the Higgs mass, and in 
particular that with least tuning (i.e., k > A/ sin 2/3), has yF so low that the NLSP should decay 
within the LHC detectors giving a displaced secondary vertex. 

The additional tuning that we discussed above for the region with k < A/ sin 2/3 is illustrated 
by the thinness of the green band in that region compared to the thicker band that we obtain for 
k > A/ sin 2/3. An interesting comment here is that the region with small k is disfavored by the 
requirement for the realistic vacuum to be a global minimum of the tree-level potential. Meanwhile 
for larger k the realistic vacuum tends to be deeper than the other minima. This matches with our 
expectation from the analytic study in Section 3.2 that larger A tends to make unrealistic minima 
deeper, while larger k may help the desired minimum to be a global one. As we can see from the 
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pink lines in the figure, in the region preferred by LHC discovery of SM-like Higgs boson, the tuning 
of the mixing quantified by eq. (19) is at the level of 1 part in 5, which is a mild one. 

The left panel in Figure 4 illustrates the situation for sizable A-terms. As in the other figures 
the green band in the plot is the region where the SM-like Higgs mass is in the range 110 - 125 GeV. 
The region of the (A, n) plane at the left of the purple dashed line is compatible with perturbative 
unification 13 . 

Apparently from the figure, there is little room to accommodate the observed Higgs boson mass 
in the region compatible with perturbative unification. Furthermore to get such large Higgs mass it 
is necessary to tune A\ in order to saturate the upper bound on the Higgs mass in eq. (56). This is 
indicated by the large values of A m j x (the pink solid lines in the figure) which shows that the region 
compatible with perturbative unification has tuning much worse than 1 part in 40. A caveat is in 
order about this result. The boundary of the perturbative region marked by purple dashed line in 
the plot should be considered as a loose one, due to the inherent uncertainties in the definition of 
a Landau pole. For instance, changing the coupling at the Landau pole by just a factor 2 would 
result in a few percent difference in the coupling A at the TeV scale. However we do not expect 
these minor corrections to change our conclusion in a significant way. 

We also remark that the region compatible with perturbative unification has a very narrow 
range of k which suggests that the phenomenology of this region of the NMSSM parameter space 
will be largely dictated by the requirement to have such a heavy Higgs boson. Note that in all the 
regions where the Higgs mass is accommodated, the realistic minimum is a global minimum. 

Finally, if we abandon the requirement of perturbative unification, then we can get the correct 
Higgs mass with much smaller fine-tuning. For example, for the choice of parameters in the figure, 
the tuning can be as low as 1 part in 10 when A > 0.7. 

6 Conclusions 

The LHC has very recently discovered what looks like a SM-like Higgs boson of mass 125 GeV. It is 
then intriguing to study implications of this observation for extensions of the SM, in particular, the 
most popular one of SUSY. It is well-known that it is difficult to accommodate such a Higgs mass in 
the minimal version of SUSY (MSSM) since the tree-level prediction for it is at most mz and loop 
corrections do not suffice to make up the difference from 125 GeV, unless we sacrifice naturalness 
of the EW scale. In this paper, we focused on a specific beyond-the-MSSM scenario, namely, the 
scale-invariant NMSSM where dimensionful terms in the superpotential can be forbidden by a Z3 
symmetry and /i-term for the Higgs doublets can be dynamically generated from the VEV of a 

13 More precisely the purple dashed line is the boundary of the region that contains couplings that stay perturbative 
(less than V2tt) up to Aqut = 1.6 • 10 16 GeV according to a one-loop RGE with boundary conditions at 1 TeV and 
4 sets of 5 + 5 of SU(5) entering in the RGE at 5 TeV. 
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singlet coupled to Higgs doublets. As a bonus of this singlet-Higgs doublet coupling, we get an 
extra quartic coupling for the Higgs doublets, which raises the tree-level SM-like Higgs mass beyond 
the mz value of the MSSM. Combined with a loop correction from stops of natural light masses it 
naively looks that the NMSSM can easily accommodate 125 GeV SM-like Higgs mass. 

However, we emphasized that there is a limitation. The above singlet-Higgs doublet coupling 
along with the VEVs inevitably results in mass terms mixing the singlet and Higgs doublet. Such 
mixing modifies the mass of SM-like Higgs from the above expectation, and the effect does not 
decouple even if the singlet is heavy. We performed a systematic study of this effect, in particular, 
we considered both the cases where the SM-like Higgs mass is reduced ("pull-down") and raised 
("push-up") as a result of the mixing. Also, we carefully distinguished between the cases with 
negligible vs. non-negligible ^4-terms and similarly, couplings being perturbative up to GUT scale 
vs. having Landau poles below. For simplicity and naturalness we chose stop mass below 500 GeV 
with small stop mixing. 

One of our sharp results is that the perturbative case which can preserve the merit of gauge 
coupling unification of the MSSM requires A\ to be tuned with the //-term in order to reach 125 GeV 
for the SM-like Higgs mass. We quantified this tuning and found that the best attainable tuning is 
roughly 20% for push-up and significantly worse tuning for pull-down. The non-perturbative option 
allows more room, for example in this case A-terms are not necessary. However we showed that 
there are still non-trivial constraints to be satisfied in order to achieve the goal of correct SM-likc 
Higgs mass. The most natural region of the non-perturbative parameter space is the one with 
k > A where tuning can be 20% and the Landau pole scale is typically below 10 4 TeV. In general, 
our findings motivate further studies of UV models which lead to desired relations among the 
parameters of the NMSSM. In particular, in the case of minimal gauge mediated SUSY breaking, 
A-terms are predicted to be small so that the perturbative case is in trouble. 

Another highlight of our work is the careful consideration of constraints from global vacuum 
stability on model parameter space. This aspect is overlooked in most of the existing literature. Our 
analysis shows this constraint can exclude significant portion of parameter space which is allowed 
otherwise. In particular, this danger of getting deeper unrealistic vacuum can be more relevant in 
the case of large A associated with low Landau pole. 

Finally we comment on some open issues which we did not elaborate in our current study. The 
current Higgs data suggests that the 125 GeV boson is roughly SM-like. We checked that in a 
sizable region of the parameter space we get a SM-like Higgs boson. In other regions of parameters 
space compatible with a 125 GeV Higgs it is possible to have modified rates for Higgs decays into 
SM states or decays into new states as neutralinos and light pseudo-scalars. However we did not 
investigate these possibilities in details. 

There are also neutral scalars and pseudo-scalars which could potentially be seen in the future 
at the LHC, although we have estimated that currently there is no significant constraint here. 
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Another example of an aspect that we did not focus on is the charged Higgs phenomenology, in 
particular, it contributes via loop along with top quark to b — >■ 57. However, this process also gets 
a contribution from chargino-stop loop, which may cancel the charged Higgs loop depending on 
parameters beyond the Higgs sector of the NMSSM, such as Wino and stop masses. Hence, we 
chose not to pursue this analysis. Finally, we did not quantify precisely the tuning of the Z mass 
due to stop contribution to rrijj as it is independent of that required to raise the tree-level mass 
of SM-like Higgs that we focused on. 

The discovery of a SM Higgs-like particle at 125 GeV has strong implications for scale-invariant 
NMSSM. We have performed a systematic analysis to explore realistic natural regions of the pa- 
rameter space where such Higgs boson can live. We find that the inhabitable space shrinks to a few 
well separated islands. This maps out well defined directions for future work both for UV model 
building and phenomenological studies. 

Note Added 

While this paper was being finalized, we learned that related ideas are being pursued independently 
by the authors of Ref. [29] . 
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